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ABSTRACT 

In this paper we present the best quahty XMM-Newton and Suzaku data from M82 X-1 so far. 
We analyze the spectra of this remarkable Ultra-Luminous X-ray Source in a self-consistent 
manner. We have disentangled emission from the host galaxy, responsible for the soft X-ray 
emission (£'<2.5keV), which is successfully described by a two-temperature thermal emis- 
sion from a hot plasma in multi-phase state, plus a narrow Fe line emission at 6.7 keV. This 
allowed us to properly study the intrinsic continuum emission from M82X-1. The continuum 
of the Suzaku spectrum is curved and the high quality data of the Suzaku spectrum allowed us 
to significantly detect a weakly broadened Fe Ka emission line. The Equivalent- Width is in 
the range 30-80 eV and it does not depend on the model applied for the continuum. Assuming 
that this line is coming from the ULX via disc fluorescence, the data indicates a disc trunca- 
tion at a radius of 6-20 gravitational radii. This value is comparable to or even larger than the 
Innermost Stable Circular Orbit of a non-spinning (Schwarzschild) black hole. Future longer 
observations might test this scenario. 
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O 1 INTRODUCTION 



, Ultra-luminous X-ray (ULX) sources are point like, ofF-nuclear ex- 
IL" • tra Galactic sources with observed luminosities greater than the Ed- 

] dington limit fo r a 10 Mq black hole (BH), with L;(>10^' erg s"' 

, ( lFabbianolll989i). The tru e nature of these objects is still open to 
j_| ■ debate jMiller & Colberlll2004h . One fundamental issue is whether 

' the emission is isotropic or beamed along our line-of-sight. A pos- 
sible scenario for geometrical beaming involves super-Eddin gton 
accre tion during phases of thermal-timescale mass transfer (King 
I2OO2I) . Alternatively, if the emission is isotropic and the Edding- 
ton limit is not violated, ULX must be fuelled by accretion onto 
Intermediate-Mass BH (IMBH), with masses in the range 100- 
10000 Mo ( Colbert & Mushotzkv 1999). Of course, it is possi- 
ble that some ULX appear very luminous due to a combination 
of moderately high mass, mild beaming and mild super-Eddington 
emission. It is also possible that ULX are an inhomogeneous pop- 
ulation, comprised of both a subsample of IMBH and moder- 
ately beamed stellar m ass black holes ( Fabbiano & White 2006 ; 
iMiller & Colberlll2004h . 

M82X-1 (also named CXOU J095550.2+694047) is one of 
the brightest ULX in the sky. The host galaxy is located nearby, 
at a distance of 3.6 Mpc (Freedma n et al.lll994 Y making M82X- 



1, with an average X-ray luminosity > 10 ergs , an excel- 
lent target for the detailed modelling of the X-ray spectrum in 
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such objects. M82X-1 is one of the few ULX for which (10- 
200 mHz ) Quasi-Pe riodic Oscillations (QPO) have been found 
Jstrohmaver & Musho tzkv 2003 ) and is th e prototype ULX with 
curved X-ray spectra. IMivawaki et al.l j200^) analyzed the Suzaku 
spectrum of this source (with an energy coverage extending up 
to to 20keV) and tested the data to a variety of models, con- 
firming that this source has a spectrum curved at high-energies. 
lMucciarellietalJ ( l2006l) analyzed the longest XMM-Newton obser- 
vations available in the archive and applied models of diffuse ther- 
mal emission from the galaxy to understand the high-energy curva- 
ture of the spectrum. The authors found that the high energy curva- 
ture could indeed be explained by the presence of a third thermal 
component having very high column densities ( A^;/~10^^ cm~^). 
This contradicts recent CHANDRA observations bv lFeng & Kaaretl 
(2010), who found column densities of Nh~IxW^^ cuT^, in fitting 
the X-ray spectra at the position of this ULX. 

As in the case of black hole binaries (BHBs), some ULXs un- 
dergo spectral transitions from a power-law dominated state to a 
"high-soft" state (see e.g. McClintock & Remillard for a descrip- 
tion of the spectral states in BHBs). During the power-law state 
the high-energy spectra of ULXs is inferred to have a power-law 
spectral shape in the 3-8 keV spectral range, together with a high- 
energy turn-over, and a soft-excess (e.g. Kaaret et al. 2006). This 
soft-energy excess can be modelled by emission coming from the 
inner-accreti on disc and is characterize d by a low disc temperature 
of !«0.2keV jMiller et al.ll2003l. l2004h . If these observed charac- 
teristics are entirely due to a BH accreting at sub-Eddington rates. 
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they imply the presence of an IMBH with M«100 - 1000 Mq. In 
the case of M82X-1, due to the huge absorption at soft X-rays by 
the host galaxy, this soft X-ray excess is not observed in the spectra. 
Syste matic studies of the hi ghest quality XMM-Newto n ULX spec- 
tra bv lYoshidaetalJ boioh and lStobbart et alj j2006h have shown 
that the majority of ULXs display a break/turnover at high energies 
(>3 k eV). Such breaks are n ot commonly seen in the spectra of 
BHBs. 'Gladsto ne et al.l ( |20 09) argue that this turnover could be due 
to Comptonization from an optically thick corona, which would 
shroud the inner regions of the accretion disc and artificially lower 
the inner temperatures obtained from simple spectral modelling. 
Under such interpretation most ULXs would then represent a new 
accretion sta te for stellar-mass BHBs. An altern ative explanation is 
proposed bv lCaballero-Garcia & Fabiad i2010l) . where the authors 
note that in many cases the turnover occurs between 5-7 keV, and 
could be due to a combination of the Iron K„ emission line and ab- 
sorption edge (6.4 and 7.1 keV respectively for neutral Iron) in a 
relativistically blurred reflection spectrum from the inner accretion 
disc of a spinning black hole. Here reflection refers to the back- 
scattering and fluorescence of X-rays dOeorge & Fabianll99lh . Un- 
der this interpretation, the soft excess emission often assumed to 
come directly from the accretion disc is also due to blurred re- 
flection, meaning that the disc emission is not directly observed. 
Other explanations have also been proposed, often based on emis- 
sion from the aforementioned slim discs in which advection of ra- 
diation from the inner regions of the accretion disc reduces the ob- 
served luminosity (Abramowicz et al. 1988). In the case of M82 X- 
1, the huge intrinsic column density of the galaxy (A'/y»10^' cm"^), 
makes the detection of soft X-ray spectral features from this source 
very unlikely. 

Here we use the best publicly available data from Suzaku and 
XMM-Newton, currently being the best available X-ray data from 
this source. Our goal is to disentangle both the emission from the 
host galaxy and the intrinsic spectrum from the source (Section[3j 
and to analyze the latter with current models used for the descrip- 
tion of the spectra from ULXs (Section|4j. We notice the presence 
of a broad Fe K line in the residuals and explore its significance in 
Section [S] Eventually, physical implications are discussed in Sec- 
tion|6] 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 XMM-Newton 

In this work, we consider the longest publicly available XMM- 
Newton EPIC-pn spectra from M82X-1 (Obs. IDs: 0112290201 
and 0206080101). The source was very bright, with a count 
rate of 3.06±0.012, 2.341+0.007 ctsjs. The datasets were obtained 
through the XMM-Newton public data archive. The EPIC-pn cam- 
era has a higher effective area than the EPIC-MOS cameras, and 
drives the results of any joint spectral analysis. The reduction and 
analysis reported in this work used SAS version 9.0.0. We checked 
for pile-up in all the observations and found that this was not sig- 
nificant (i.e. it is <5% for the high energy channels) for both obser- 
vations. 

In Table [T] we present a log of the observations. We applied 
the standard time and flare filtering (rejecting high-background pe- 
riods of rate >0.4 counts/s, as recommended for the pn camera[j). 



Information provided at "node52.html" of the User Scientific Guide. 



We filtered the event files, selecting only the best-calibrated events 
(pattem<4 for the pn), and rejecting flagged events (flag= 0). 

For each exposure, we extracted the flux from a circular region 
on the source (centre at coordinates RA = 9 /i55 ra50.2 s deg, Dec = 
69 (jf40 ra47 .s deg and radius 18arcsec). The background was ex- 
tracted from a circular region, not far from the source and away 
from boundaries of the chips and the nucleus of the galaxy. We 
built response functions with the SAS tasks rmfgen and arfgen. 
We fitted the back ground-subtrac ted spectra with standard models 
in XSPEC 12.5.0 ( lAmaudllT99^ . All errors quoted in this work 
are 90% confidence errors. For the spectral fitting we used the 0.5- 
lOkeV range. The resulting spectra were grouped with the FTOOL 
grppha to bins with a minimum of 20 counts each. 

2.2 Suzaku 

M82X-1 was observed with SUZAKU jMitsuda et all l2007h on 
three occasions in 2005 October 4, 19 and 27 for approximately 
32, 40 and 28 ks respectively. The four operating detectors consti- 
tuting the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
were operated in the "normal" clock mode in all observations. In 
all cases, the detectors were operated in both the 3x3 and 5x5 
editing modes. A total (co-added) front-illuminated (FI) exposure 
of approximately 96, 110 and 80 ks were obtained for the three 
respective observations. The corresponding back-illuminated (BI) 
exposures were approximately 32, 37 and 27 ks. Using the lat- 
est HEASOFT software package we processed the unfiltered event 
files for each of the XIS CCDs and editing modes operational in 
the respective observations, following the SUZAKU Data Reduc- 
tion Guid^ We started by creating new cleaned event files by re- 
running the SUZAKU pipeline with the latest calibration, as well 
as the associated screening criteria files. XSELECT was used to ex- 
tract spectral products from these event files. In all observations, 
source events were extracted from a circular region of 180" radius 
centred on the point source, and background spectra from another 
region of the same size, devoid of any obvious contaminating emis- 
sion, elsewhere on the same chip. The script "xisresp'0 with the 
"medium" input was used to obtain individual ancillary response 
files (arfs) and redistribution matrix files (rmfs). "xisresp" calls 
the tools "xisrmfgen" and "xissimarfgen". Finally, we combined 
the spectra and response files from the three front-illuminated in- 
struments (XISO, 2 and XIS3) using the FTOOL ADDASCASPEC . 
This procedure was repeated for each observation resulting in a 
total of six XIS spectra. After inspecting that there were no signifi- 
cant spectral differences between the three observations we further 
co-added the various front and back-illuminated spectra. This re- 
sulted in a total good exposure of approximately 280 and 92 ks for 
the front and back-illuminated spectrum respectively. Finally, the 
FTOOL GRPPHA was used to give at least 20 counts per spectral 
bin. The two XIS spectra was fit in the 1. 5-10.0 keV energy range 
with the 1.7-2.1 keV energy range being ignored due to the possi- 
ble presence of un-modeled instrumental features^ 

For the hard X-ray detector (HXD. iTakahashi et al.ll2007h we 
again reprocessed the unfiltered event files for the respective ob- 
servations following the data reduction guide (only the PIN data is 
used in this analyses). Since the HXD is a collimating rather than 
an imaging instrument, estimating the background requires individ- 
ual consideration of the non X-ray instrumental background (NXB) 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/ 

^ http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/xisresp 
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Table 1. XMM-Newton and Suzaku observations log. 



Satellite 


Obs 


Obs ID 


Date 


Exposure time ( s) 


XMM-Newton 


1 


0112290201 


2001-05-06 


30558 


XMM-Newton 


2 


0206080101 


2004-04-21 


104 353 


Suzaku 


3 


100033010 


2005-10-04 


32 327 


Suzaku 


3 


100033020 


2005-10-19 


40358 


Suzaku 


3 


100033030 


2005-10-27 


28 363 



and cosmic X-ray background (CXB). The appropriate response 
and NXB files were downloaded for the respective observation^; 
in each case the tuned (Model D) background was used. Common 
good time intervals were obtained with MGTIME which combines 
the good times of the event and background files, and XSELECT 
was used to extract spectral products. Dead time corrections were 
applied with HXDDTCOR, and the exposures of the NXB spectra 
were increased by a factor of ten, as instructed by the data reduc- 
tion guide. The contribution from the CXB was simulated using 
the form of Boldt ( 1987), with the appropriate normalisation for the 
nominal pointing (all observations were performed with HXD nom- 
inal pointing), resulting in a CXB rate of ^ 2.97 x 10"^ count s"'. 
The NXB and CXB spectra were then combined using MATHPHA 
to give a total background spectrum, to which a 2% systematic un- 
certainty was added. Similarly to the XIS data, we co-added all 
three PIN observations and grouped the final spectrum to have 
a minimum of 500 counts per energy bin to improve statistics, 
and again allow the use of minimization during spectral fitting. 
The PIN data reduction yielded a total (co-added) source rate of 
(1.5 ± 0.3) X 10"- counts"' which is 2.5% of the total observed 
flux, with a good exposure time of « 93 ks. 



3 THE DETERMINATION OF THE EMISSION FROM 
THE DIFFUSE COMPONENT 

The soft energies in the spectra of M82X-1 are dominated by 
the presence of X-ray emission lines having constant energies and 
width, thus indicating an origin in the diffuse gas of the galaxy, 
rather than intrinsic to M82X-1. We performed simultaneous fits 
of all the observations, considering constant properties (in time) for 
both thermal components of the diffuse emission from the galaxy 
(i.e. temperatures and element abundances tied between the obser- 
vations, although see below for the latter). This is as expected, since 
the physical properties of the galaxy are likely to be constant with 
time. To account for these lines we initially applied a single ther- 
mal plasma model with variable metal abundances (vmekal model 
in XSPEC) to describe the softest part of the spectrum (< 2.5 keV). 
However, strong residuals were left with this model and we thus 
added a further second thermal component (as previously done in 
iMucciarelli et al. |[2003) resulting in a significant improvement to 
the fit, with the residuals at < 2.5 keV being almost null, within 
the instrumental errors. The values obtained for the spectral pa- 
rameters are shown in Table [2] The valu e for the temperature of 
the hottest component in Obs. 2 found by IMucciarelli et alj bOOd) 
{kT - 1 .67+0. 1 1 keV), is different possibly due to the addition of a 
further thermal disc component in their analyses. We tested for the 
presence of this further component and did not find it to be signifi- 
cant in the spectra. 

http://www.astro.isas.ac.jp/suzaku/analysis/hxd/ 



Table 2. Best-fit model spectral parameters with 90% confidence errors ob- 
tained for the Global-Simultaneous fit shown in Figure[T] 



Models 


Parameters 


Obs 1 


Obs 2 


Obs 3 


Diffuse 




1.42±0.03 


1.25±0.03 


0.58±0.01 


Plasma 


(10^2 cm-2) 










kT(l)(keV) 


= Obs3 


= Obs3 


0.96±0.02 






0.13±0.02 


0.18±0.02 


0.06±0.02 




(10^2 cm-2) 










kT(2) (keV) 


= Obs3 


= Obs3 


0.60±0.02 


Narrow Fe 


E| (keV) 


= Obs3 


= Obs3 


6.7±0.3 


Line 


o- (keV) 


= Obs3 


= Obs3 


0.06±0.04 


Highecut 


r 


1.47±0.03 


1.49±0.03 


1.79±0.02 


Powerlaw 


Ec (keV) 






6.83±0.19 




Ef (keV) 






9.7±1.4 



We fixed metal abundances to solar values for Obs. 1 and 
2 (XMM-Newto n spectra), in agreeme nt with previous work on 
XMM/RGS data teead & Stevendl 20021) . Setting solar metal abun- 
dances for Obs. 3 (Suzaku spectrum) provided a very poor fit 
(X^/v = 1.92, V = 3450) and visible low energy residuals. The 
best fit was achieved allowing metal abundances free to vary in 
the Suzaku spectra and they turned out to be super-solar for certain 
elements (Z = 6 - BxZ© for Si, S and Fe). This fit apparently indi- 
cates super-solar abundances for the diffuse emission and therefore 
it is ha rd to reconcile with previous studies o f metallicity of this 
galaxy iOriglia et all2004lRanalli et alj|2008l) . These lines might 
have a difi^erent origin apart from the emission from the relatively 
cold emission from the galaxy (nevertheless, we refer to this emis- 
sion as hot hereafter), which the model can not account for. Further 
studies of these emission lines require the use of additional emis- 
sion processes/components that we will not consider in the present 
paper. 

There is also a narrow line (i.e. cr = 0) centred at a;6.7keV 
coincident in energy with Fe He„. We included it in all the fits re- 
ported above and in the following fits as well. This line can not be 
accounted for the two thermal models (with solar abundances) em- 
ployed to model the soft X-ray emission. We extracted spectra from 
other regions of the galaxy and confirm the presence of the narrow 
line. Stri ckland & Hec kman ( 2007) spectrally resolved the Fe line 
at 6.7 keV with Chandra and XMM-Newton data and showed that 
it is primarily distributed diffusely rather than associated with dis- 
crete sources. 

Having confirmed that the combination of two thermal plasma 
components plus emission from a narrow Fe line indeed provides 
a first approximation description of the emission from the diffuse 
component in the global-simultaneous fits, we proceed with this 
combination in all fits presented hereafter. As previously noticed by 
iMivawaki et al.l 1(2009) for Obs. 3 (for which the energy coverage is 
maximum 0.5 - 20 keV in contrast to 0.5 - 10 keV for the remain- 
der observations), the high-energy spectrum is strongly curved. We 
applied a curved phenomenological model to describe the high- 
energy spectrum (high-energy cut-off in XSPEC). From the 
global simultaneous fits , the cut-off for Obs . 3 was constrained to 
be Ec = 6.83±0.19keV. lMizuno et al . I ( |2007|) reported a very simi- 
lar cut-off at 6-7 keV in the spectrum of NGC 1 3 1 3 X-1 . The results 
of the global-simultaneous fits are shown in Figure[T]and Table|2l 
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M82 X-1 




Energy (keV) 

Figure 1. Global-Simultaneous fit of spectra of Obs. 1, 2 and 3 with a 
continuum phenomenological model (Highecut Powerlaw for the high- 
energy emission). Dift'use emission from two hot plasmas (kT = 0.6 keV 
and kT = 0.96 keV) and solar metallicities for XMM-Newton spectra and 
emission from a narrow Fe line have been taken into account. Black, 
red, blue and green correspond to FI,BI Suzaku/XIS, XMM-Newton and 
Suzaku/PIN/HXD spectra, respectively. 

4 THE DETERMINATION OF THE BEST MODEL FOR 
THE CONTINUUM INTRINSIC TO THE SOURCE 

Considering the results obtained in the previous Section regarding 
the X-ray emission from the diffuse emission of the galaxy, we are 
now in good condition to determine the best model for the intrinsic 
emission of the source. We can see in Figure [T] and in the power- 
law photon indices obtained (see Table O that the intrinsic spec- 
trum of M82X-lis variable. Thus we fitted each spectrum sepa- 
rately and we froze the parameters regarding the emission from the 
diffuse component from the model of the previous Section (i.e. 2 
mekal components plus narrow Fe line), expected to be constant. 
To this base-line model we added further models to describe the 
high-energy emission, as explained in the following. 

The high-energy spectra are strongly curved, and we find 
that a model consisting by an absorbed power-law is not a good 
description of the data of Obs. 3 (y^/v = 1.25, 1.28,2.16 with 
V = 1340, 1399, 730 for Obs. 1, 2 and 3, respectively). The results 
from the Suzaku spect rum are in agreement with the analysis from 
iMivawaki et al. Il l2009l) . We considered X-ray emission from a mul- 
ticolor disc (diskbb in XSPEC) and Comptonization (CompTT in 
XPSEC). 

The spectra of M82X-1 are con vex, representing very well 
the prototype of ULXs of this class ( IMakishimal l2007h . Indeed, 
the spectra are well fitted with emission from a rather hot inner 
disc (kTi„ = 2.7 - 3.6 keV; see Table[3ll, values in agreement with 
previ ous studies jStrohmaver & Mushotzkvll2003l : lMivawaki et al.l 
I2OO9I) . Nevertheless, these values are significantly hotter than those 
obtained in recent ob servation (kTi„ = 0.9 - 1.6) of M82X-1 by 
iFeng & Kaare3 l l2O10l) . This fact might correspond to the known 

i iositive temperature-lu minosity correlation known to act in ULXs 
Makishima et al .I2OOOI) . But we do not observe a tr ue change in the 
lumino sities of M82X-1 with respect those from iFeng & Kaared 
(l2010l). being very s imilar to the ones reported in this paper. 
iFeng & Kaaretl ( I2OIO I) identify their spectra as being from M82 X- 
1 in a state similar to the canonical high-soft state in BHBs (i.e. 
they follow the relationship LocT*). M82X-1 might indeed be in 
a diffe rent state during these observations, instead. iMizuno et al.l 
( I2OO7I) reported the cut-off in the spectrum of NGC 1313 X-1 in 
the power-law state of ULXs and IMivawaki et"ZI ( l2009b argued 



Table 3. Models applied for the continuum intrinsic to M82X-l(i.e. with- 
out contribution from the galaxy). Intrinsic luminosities are calculated in 
the 0.5-10 keV and 2-20 keV energy range for XMM-Newton and Suzaku 
spectra, respectively. 



Models 


Parameters 


Obs 1 


Obs 2 


Obs 3 


Diskbb 


kTi„ (keV) 


3.49±0.16 


3.05±0.07 


2.73±0.06 






1.19 


1.17 


0.82 






(1589/1340) 


(1636/1399) 


(465/566) 




LxCxlO*") 


3.22±0.02 


2.05±0.20 


4.8±0.4 




(ergs"') 








CompTT 


kTe (keV) 


2.08±0.02 


2.21±0.02 


2.26±0.02 




T 


22.0±1.1 


10.1±0.2 


8.1±0.3 






1.05 


1.14 


0.82 






(1400/1338) 


(1589/1397) 


(465/564) 




LxCxlO'*") 


6.6±3.0 


4.15±2.0 


4.7±0.4 



(ergs ') 



that M82X-1 was in the power-law state during Obs. 3, as well. 
Since there is barely any spectral evolution during these observa- 
tions, we conclude that the source was in the power-law state during 
all the observations reported in this paper. 

The X-ray spectra of M82X-1 are reproduced very well by 
the compTT model (see results of spectral fits in Table [3]l, invok- 
ing a rather low electron temperature (kTe = 2.0 - 2.3 keV) and a 
large optical depth (t = 7 - 23). The results of Obs. 3 agree with 
those obtained previously for M82X-1. The statistics reported in 
this work for Obs. 3 is poor ix^ lv>\.2), since we considered the 
full 0.5-20 keV energy range in the fits of the Suzaku spectrum and 
below 2 keV there are unavoidable residuals most probably due to 
calibration effects (and/or) the limitation of our model to describe 
the diffuse emission from the galaxy. Since our goal is the study of 
the high-energy emission intrinsic to the source we will just con- 
sider the 2-20 keV energy range in the fits reported hereafter. 

We can see from the results of fitting with both X-ray curved 
models that the description with a Comptonization model is a bet- 
ter description just for the spectrum of Obs. 1 (significant improve- 
ment by 8cr). We thus consider the Comptonization model a better 
description of the spectrum of Obs. 1. If luminosities from M82 X- 
1 are truly super-Eddington, then the m entioned paramet ers suggest 
that the cold-corona model proposed bv lGladstone et al] | [200^) is a 
possibility in the description of the spectra of ULXs, emitting in the 
so-called Very-High state. Nevertheless, this statement is just tenta- 
tive and it is not required for the full data set analyzed in this paper. 
We will then consider both curved models (diskbb and compTT) 
as possible descriptions of the data and we will not discuss their 
validity any further. 



5 CONSTRAINTS OF AN EMISSION BROAD FE Ka 
LINE IN THE SUZAKU SPECTRUM 

The results of the previous Section with feasible models for the 
continuum leave some positive residuals in the 6-7 keV energy 
range in the SuzakufXlS spectrum of Obs. 3 (see left panel of Fig- 
ure [3}. These residuals could not be accounted for by the addition 
of the narrow Fe line and indicate the likely presence of a broad 
Fe emission feature. To analyze this possibility further, we added a 
broad gaussian line to the best fit model (described in the previous 
Sections) of Obs. 3. 

The addition of a gaussian centred at 6.4-6.97 keV in the best 
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M82 X-1 (Obs. 1) 




EriGrgy (keV) 
M82 X-1 (Obs. 2) 





IlliiiljIiiiilJiiii 



Energy (keV) 
M82 X-1 (Obs. 3) 



I'll Miiii44* 




■f Till I 



EnGrgy (keV) 

Figure 2. XMM-Newton and Suzaku spectra of M82X-1 (upper to lower) 
fitted with compTT model. Continuous curvature of the overall emission and 
diffuse emission from two hot plasmas (kT = 0.6 keV and kT = 0.96 keV) 
and galactic emission in the form of a narrow Fe line have been taken into 
account. 



model describing the spectrum of Obs. 3 improves significantly the 
quality of the fit with any of the models used to describe the hard 
02 kev) X-ray spectra (improvement of Ax^ = 13, 22 for v = 3 ad- 
ditional degrees of freedom for Disc and Comptonization models, 
respectively). This means a 3.3 - 4.4cr detection, which means a F- 
test probability of < 0.1% that the improvement occurs by chance 
with any of both models. 

The line is broad (tr = 0.3 - 0.4 keV) and the derived equiva- 
lent width (EW) of the line is of EW 



5St_ll, 59t_lj eV for the Disc 



and the Comptonization models, respectively. We conclude that the 
value of the equivalent width of the line is significantly different 
from zero and that it is model-independent. Thus, we can consider 
the detection of the broad Fe line as a robust one, which does not 
depend upon the model considered for the continuum. In the right 
panel of Figure |3] we plot both the broad Fe line and the best de- 
scription used for the underlying continuum. 



5.1 Constraints on the spin of the Black Hole 

As shown above the presence of a broad Fe line is significant in 
the Suzaku spectrum (Obs. 3). The Equivalent-Width is in the range 
30-80 eV, in agreement with findings for some Active Galactic Nu- 
clei, and it does not depend on the model applied for the contin- 
uum (Disc emission or Comptonization models). Accretion mod- 
els predict that fluorescence lines broadened by relativistic effects 
could arise from reflection of X-ray emission onto the inner region 
of the accretion disc surrounding the black hole. In order to de- 
rive the inner disc radius and inclination we fitted the spectrum of 
Obs. 3 with a relativistic disc line profile calculated in a maxim ally 
spinning Kerr space-time (laor model in XS PEC : iLaoill 1 99 ih but 
with the inner radius allowed to vary. We froze the emissivity in- 
dex of the power-law emission and the outer radius to q = 3 and 
Rout = 100 Rg, thus canonical values in BH-accreting systems. The 
inner disc radius was found to be Ri„ = 1 1+^, lO^J Rg (values corre- 
sponding to the Disc and Comptonization models, respectively, for 
the continuum). This value is comparable to or even larger than the 
innermost stable circular orbit of a non-rotating (Schwarzschild) 
BH (i.e. 6Rg). The inner disc inclination was found to be very 
small, too (i = IS^'j, 10^]g°), thus implying a close face-on con- 
figuration. Adding this line supposed an improvement in the fit of 
Ax' = 15,21 for V = 4 additional degrees of freedom for Disc and 
Comptonization models, respectively, thus very similar to the case 
of the inclusion of the broad gaussian. 



6 DISCUSSION AND CONCLUSIONS 

In the present work we have studied the best broad-band X- 
ray spectra of M82 X-1 currently available in X-rays from XMM- 
Newton and Suzaku satellites. With these high-quality data we have 
been able to disentangle emission from the diffuse plasma of the 
galaxy and intrinsic emission from the ULX. Our results imply 
that the X-ray diffuse emission is composed by a two-temperature 
thermal emission from a hot plasma in multi-phase state. Similar 
temperatures of the components {kT = 0.96±0.02, 0.60±0.02 keV) 
have been fou nd for the highest S/N ASCA spectra of nearby ellip- 
tical galaxies jBuote & Fabianll998h . The intrinsic spectrum of the 
source obtained with Suzaku is significantly curved and properly 
accounted by models with intrinsic curvature (either emission from 
the inner accretion disc or Comptonization from cold electrons in 
an optically thick corona), in agreement with previous studies. But 
the most important result of this paper is that the (very) high qual- 
ity observation with Suzaku has allowed us to detect a significant 
broad Fe Ka line in the spectrum of this notorious ULX. Assuming 
that the line feature is solely coming from the ULX via disc fluo- 
rescence, and that the inner radius reaches the last stable orbit, the 
data rules out a rapidly spinning BH and favours a non-sp inning 
Schwarzschild-BH. Previous studies ( lFeng&Ka^l2010i) have 
suggested that the BH in this ULX is rapidly rotating but the possi- 
bility of a non-spinning BH has not been rejected so far. 
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The obtained equivalent width (EW) is in the range SO- 
SO eV. This value for the EW (in Obs. 3) is notoriously lower 
than the t entative value previously reported (EW=0. 26-0.43 keV; 
IStrohmay er & Mushotzky 2003) for the spectrum of Obs. 1. With 
the application of our model for the diffuse emission and the contin- 
uum of the source, we do not detect any Fe line in the same dataset 
of their work (Obs. 1). We have found that in the Suzaku spectrum 
this line is broad and the EW is similar to typical values found for 
the b road Fe K line in a wide sample of AGN jde la Calle et al.l 
The value we found for the EW is significantly smaller than 
those found in the typical c ases of claime d relativistic lines (e.g., 
«300 eV in MCG-6-30-15; iMiniutti et ai . 2007), for which an ex- 
tremely fast rotating Kerr-BH has been clai med. Nevertheless , re- 
cent studies are setting much lower values jNoda et al]|201ll) for 
the spin of this AGN, underlying in the fact that the continuum was 
poorly modelled. 

In the previous w ork of analysis of the Suzaku spectrum by 
iMivawaki et al. the broad Fe line was described by the ad- 

dition of several narrow lines. This hypothesis looks reasonable, 
sin ce a similar complex has been clearly seen in our Galactic Cen- 
tre jKovama et akluOOTh . Nevertheless, the presence of relativisti- 
cally broad Fe lines seems to be common in all types of accreting 
compact objects (neutron stars, BHB, AGN) and evidence of their 
finding in ULXs would be very helpful in the understanding of the 
real physical processes accounting in these sources. Future longer 
X-ray observations will certainly provide further insights into the 
Fe complex and the presence of additional reflection features for 
this archetype of the ULX class. 
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